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Quantum Heisenberg antiferromagnets in a uniform magnetic field: Non-analytic 
magnetic field dependence of the magnon spectrum 
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We re-examine the l/S'-correction to the self-energy of the gapless magnon of a D-dimensional 
quantum Heisenberg antiferromagnet in a uniform magnetic field h using a hybrid approach between 
1/S-expansion and non-linear sigma model, where the Holstein-Primakoff bosons are expressed in 
terms of Hermitian field operators representing the uniform and the staggered components of the 
spin-operators [N. Hasselmann and P. Kopietz, Europhys. Lett. 74, 1067 (2006)]. By integrating 
over the field associated with the uniform spin-fluctuations we obtain the effective action for the 
staggered spin-fluctuations on the lattice, which contains fluctuations on all length scales and does 
not have the cutoff ambiguities of the non-linear sigma model. We show that in dimensions D < 3 the 
magnetic field dependence of the spin- wave velocity c_ (h) is non-analytic in h 2 , with c_ (h) — c— (0) oc 
h 2 In \h\ in D — 3, and C-(h) — c_(0) oc \h\ in D — 2. The frequency dependent magnon self-energy 
is found to exhibit an even more singular magnetic field dependence, implying a strong momentum 
dependence of the quasi-particle residue of the gapless magnon. We also discuss the problem of 
spontaneous magnon decay and show that in D > 1 dimensions the damping of magnons with 



momentum k is proportional to ft 



if spontaneous magnon decay is kinematically allowed. 
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I. INTRODUCTION 



One of the most successful methods for obtaining the 
low-temperature properties of ordered quantum Heisen- 
berg magnets is the expansion in inverse powers of the 
spin quantum number S. The idea is to first map the 
spin Hamiltonian onto an interacting boson model us- 
ing either the Holstcin-PrimakofS- or the Dyson-Maleyev 
transformation 2 '-, and then study the resulting inter- 
acting boson system by means of the usual many-body 
machinery. As the interaction vertices appearing in the 
boson Hamiltonian involve the small parameter of 1/S, 
the perturbative treatment of the interaction is formally 
justified for large S. See, for example, Refs. [jQ for 
early applications of this approach to quantum antifer- 
romagnets (QAFM). A disadvantage of this method is 
that calculations for QAFM beyond the leading order in 
1/5 are very tedious due to a large number of interac- 
tion vertices^. Moreover, the vertices are even singular 
for certain combinations of external momenta££i Al- 
though the singularities cancel in physical quantities if 
the total spin is conserved^, the appearance of singular- 
ities at intermediate stages of the calculation indicates 
that this approach is not always the best way of calcu- 
lating fluctuation corrections to the magnon spectrum. 

In this work we shall re-consider the leading l/S'- 
correction to the magnon self-energy of spin-S quantum 
Heisenberg antiferromagnets in a uniform magnetic field 
h at zero temperature in the regime where the system 
has a finite staggered magnetization. Our starting point 



is the Heisenberg Hamiltonian 



H 



2_^ Jij Si ■ Sj 2_^ " ' °i 



(1.1) 



where Si are spin operators normalized such that Sf = 
S(S + 1) and the magnetic field h is measured in units 
of energy. The exchange integrals Jy connect nearest 
neighbor sites r, and Tj on a D-dimensional hypercubic 
lattice with lattice spacing a, total volume V = a D N and 
N sites. As long as \h\ is smaller than a certain critical 
value h c (see Eq. (|2.20[) below), the spin configuration 
in the ground state is canted, as shown in Fig. [TJ We 
choose our coordinate system such that the magnetic field 
h = he x points along the x-axis and the staggered mag- 
netization M s — M s e z points in z-direction. The mag- 
netic field generates a uniform magnetization M = Me x 
pointing in the same direction as h, giving via h a gap in 
the transverse magnon polarized parallel to h, while the 
magnon polarized perpendicular to h remains gapless. 

Due to the canting of the spins, the effective boson 
Hamiltonian obtained from Eq. (|1.1|) within the Holstein- 
Primakoff transformation contains cubic interaction ver- 
tices proportional to S 1-1 / 2 . Hence, to obtain the com- 
plete l/S'-correction to physical observables, the cubic 
vertices should be treated in second order perturbation 
theory. The leading 1/S-corrections to the magnon spec- 
trum turns out to be rather peculiar: Zhitomirsky and 
Chcrnyshev^ have shown that for intermediate magnetic 
fields in a certain range h* < \h\ < h c there are no well- 
defined magnons in a large part of the Brillouin zone due 
to spontaneous two-magnon decays. Moreover, Syromy- 
atnikov and Maleyev-ii. calculated the l/S'-correction to 




m 



FIG. 1: (Color online) Spin configuration (Si) = Srhi in the 
classical ground state of a two-sublattice antiferromagnet sub- 
ject to a uniform magnetic field h = he x in the ^-direction. 
The hypercubic lattice can be divided into two sublattices, la- 
beled A and B, such that the nearest neighbors of a given site 
all belong to the other sublattice. The solid square denotes a 
site of the A-sublattice and a solid circle denotes a site of the 
B-sublattice. Here #o is the classical canting angle between 
the direction of the staggered magnetization fli and the local 
spin direction Aj. 



the anisotropy induced gap of the magnon polarized par- 
allel to the magnetic field, and showed that in dimensions 
D < 3 the correction is unexpectedly large. They sug- 
gested that meaningful results can only be obtained if 
the 1/S-expansion is re-summed to all orders, which is of 
course impossible in practice. 

Unfortunately, within the conventional 1/S expan- 
sion, the expressions for the magnon self-energies (see 
Refs. [lOlflll ]) are quite complicated. For example, from 
the expression for the magnon self-energy given by Zhit- 
omirsky and Chernyshev 10 (which we reproduce in Ap- 
pendix B) it is not immediately obvious that one of the 
magnon branches remains gapless. In this work we shall 
therefore re-consider this problem using our recently pro- 
posed parameterization of the l/S'-expansion in terms 
of Hermitian field operators^. The advantages of such 
an approach have already been pointed out in Ref. |7j, 
but the practical usefulness of this method has not been 
demonstrated. In a sense, our method is a hybrid ap- 
proach between the l/S'-expansion and the non-linear 
sigma model (NLSM) approac h 9 i 12 i 13 . Recall that the 
NLSM is an effective continuum theory for the staggered 
spin-fluctuations of a QAFM. In contrast to the singu- 



lar interaction vertices encountered in the conventional 
1/S-expansion, the vertices describing interactions be- 
tween transverse spin-fluctuations in the NLSM are finite 
in momentum space and all scale as k 2 for h = 0. On 
the other hand, the NLSM has to be regularized using an 
ultraviolet cutoff, so that the NLSM approach cannot be 
used to obtain the numerical value of observables which 
receive contributions from wave- vectors in the entire Bril- 
louin zone. Our approach combines the advantages of the 
1/S-expansion with the those of the NLSM by parame- 
terizing the degrees of freedom in the 1/S-expansion from 
the beginning in terms of a lattice version of the contin- 
uum field representing staggered spin fluctuations in the 
NLSM. 

The rest of this work is organized as follows: After 
giving a detailed description of our hybrid approach in 
Sec. HH we derive the effective action for staggered spin 
fluctuations of our lattice model in Sec. IIIII and exhibit 
the precise connection with the NLSM, where only the 
leading orders in the derivatives are retained. In par- 
ticular, we show how the regular vertices of the NLSM 
emerge from the conventional 1/S-expansion. In Sec. IIVI 
we then use our method to derive expressions for the fre- 
quency dependent part of the magnon self-energies which 
for small magnetic field h determines the dominant fa- 
dependence of the magnon dispersions. In Scc.|V]thc self- 
energy of the gapless magnon is evaluated; in particular, 
we show that in dimensions D < 3 the fluctuation cor- 
rections to the spin-wave velocity and the quasi-particlc 
residue of the gapless magnon exhibit a non-analytic h- 
dependence. We also discuss the problem of spontaneous 
magnon decay in general dimensions. After a brief sum- 
mary of our results in Sec. IVI1 we give in Appendix A 
explicit expressions for the quartic interaction vertices 
associated with two-magnon scattering in our hybrid ap- 
proach. Finally in Appendix B we show numerically that 
in D = 2 our result for the magnetic field dependency of 
the spin-wave velocity of the gapless magnon can also be 
extracted from the self-energy given by Zhitomirsky and 
Chernyshev in Ref. [l(| . 



II. HYBRID APPROACH: COMBINING THE 

ADVANTAGES OF THE 1/S-EXPANSION WITH 

THOSE OF THE NLSM 

A. Holstein-Primakoff boson Hamiltonian 

For completeness, let us briefly recall the general pro- 
cedure for setting up the 1/S-expansion around a given 
classical ground-state, characterized by the directions 
rhi = (Si)/\(Si)\ of the local magnetic moments^. Sup- 
plementing the unit vector rhi by two additional unit 
vectors e\ ' and e\ ' such that e\ ,e\ , rhi form a right- 
handed orthogonal triad of unit vectors, and defining the 
corresponding spherical basis vectors eP i = e\ + ipe\ , 
p = ±, we express the components of the spin operator 
Si in terms of canonical boson operators hi and b\ using 



the Holstein-Primakoff transformation^, 

5, = Sfrh t + S± = s}rh t + ~ £ S^ 



p c p 



D = ± 



with 



S] — 5 - rii , n 4 = b\bi , 



2S l ~W' 



256J Wl 
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(2.1) 

(2.2a) 
(2.2b) 

(2.2c) 



Our spin Hamiltonian (jl.lj) can then be written as the 
following bosonic many-body Hamiltonian^ 



H = Ef + H\ +H l }+H ± + H' 



with the classical ground state energy 

S 2 



and 



E ° " -J- Yl J ' ; : 



s 



-Sj> 



m,i, 



(2.3) 



(2.4) 



h\ = -g- ^~1 Jijihi ■ rhj (nj + nj) + y^ y h- m^ , (2.5) 



H = 



H- 






(2.6) 



(2.7) 



y pp' 



-i"*ij 



» 3 

= -YjJij{Si--m j )n j 

ij 

-J2s^{h-Y,JvSrh 3 ). (2.8) 



The part H' of the Hamiltonian describes the coupling 
between transverse and longitudinal spin fluctuations 
generated by the uniform magnetic field. Within the 
Holstein-Primakoff approach, we expand the square roots 
in Eqs. (|2.2bp and l[2"T2"c|) in powers of S^ 1 , 



St = V2S 



Si 



'2S 



Is" 

b\rn 

45 



(2.9a) 
(2.9b) 



The boson representation of the operator H can then 
be written as an infinite series of multiple-boson interac- 
tions involving even powers of boson operators, while H' 
becomes an infinite series of terms involving odd powers 
of boson operators, 



H- 1 = H^ + Hl + 0{S~ 
H' = ^1+^ + 0(5 



-l/2x 



(2.10) 
(2.11) 



where the subscripts indicate the number of boson oper- 
ators. Making the reasonable assumption that the true 
spin configuration in the ground state resembles the clas- 
sical one shown in Fig. [T] (but with a renormalized canting 
angle #), we have 



(Si) = \(Si)\rhi , rh, = Qne z + me x 



(2.12) 



where we have chosen h = he x , and the true canting 
angle $ is related to n and m via n = cos i? and m = sin $. 
Here Q assumes the value +1 on one sublattice (which we 
call the A-sublattice) and — 1 on the other sublattice (the 
B-sublattice). A convenient choice of the other members 
of the local triad is 



,(!) 



D (2) 



CjiTb&x 



(2.13) 



The relevant scalar products in this basis are for nearest 
neighbor sites i and j, 



(2.14a) 
(2.14b) 
(2.14c) 

-i\d, (2.14d) 
(2.14e) 
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= 2n 2 , 


et 
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= e~ 


<- 


= 2m 2 , 


et- 


rhj 


= -el 


' ■ m 


j = —2in 
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rhi 


= hm 






where we have defined 






a = 


-- n 2 - 


m = 


= 1 - 2m 2 




A = 


= 2nm 


= sir 


(21?)- 



cos(2tf), 



Then we obtain from Eq. (|2.4j) . 
E$ = -NDJS 2 a 
and from Eq. (|23| . 



NShm, 



H l 2 = -irhc^m, 



where 



Z h = l 



2mdh 



and we have introduced the notation 



h c 
Sh 



ADJS, 

h — h r m. 



(2.15) 
(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 

(2.21) 



In the classical limit S — ► oo the exchange field h c m 
exactly cancels the external field h, so that in this limit 
Sh = 0. However, for finite S the difference Sh = h — h c m 
is finite. We shall show in Sec. IHII that Sh is actually of 

the order of mh c /S. The longitudinal part H\ of the 
Hamiltonian involving four boson operators is 



m 



■£'« 



(2.22) 



Let us emphasize that if we use the Dyson-Maleyev 
transformation^ to bosonize the spin operators, we ob- 
tain a non-Hermitian transverse part H± which differs 

from Eq. (|2.24|) while Hi, H2, H3, and H\ are the same 
as above. Since the physical quantities calculated in this 
work are essentially determined by H3 our results do not 
depend on whether we use the Holstein-Primakoff or the 
Dyson-Maleyev formalism. 



and the leading two terms of the transverse part of the 
Hamiltonian are 

ij 

+ (e+ ■ e+) b\b] + (e- ■ ej) b^} 
= I £ Jio Vn 2 (b\b 3 + &}&,-) + n a (b\b] + bjh)] , 



(2.23) 



- 1 n x— * 



w. 



>,t„.fct 



n z bibj + biUjbj + b\b^nj + b]mbj 



— 5~£ Jij n ibib] + bibjiij + b\rnbj + b\n,jbj . 

ij 

(2.24) 

Finally, the part H' of our effective boson Hamiltonian 
describing the coupling between transverse and longitu- 
dinal fluctuations can be written as 

H' = \Y,JvQSl 2) n 3 +n5hJ2QSl 2) , 



where we have set Sir = S^ ± iS^ , so that 



sP - (1) 



er'-S t = -(S+ + Sr) 



s 



(2) = e (2) 



' °i — 2j(^i °i )• 



(2.25) 

(2.26a) 
(2.26b) 



The alternating factor Q in Eq. (|2.25[) indicates that this 
term describes Umklapp scattering across the boundary 
of the antiferromagnetic Brillouin zone. For our purpose 
it is sufficient to neglect all terms in the expansion of 
Eq. (|2.If [) involving five and more boson operators, which 
amounts to retaining only Hi and H3. With our choice 
of basis vectors these can be written as 



(2.27) 



6 1 =nSh^— 1 £ i Ci(bi-bt), 



H, 



2~ £-JyCi(&»-&J>j 



(2.28) 



B. Linear spin- wave theory 

To obtain the magnon spectrum within linear spin- 
wave theory, we neglect HJ and H', and approximate the 
transverse part H by its quadratic term in the expan- 
sion of the spin operators in terms of the boson operators, 
H^- s=s H2 ■ We should now diagonalize the quadratic 

boson Hamiltonian Hi = H\ + H^- We work in the 
sublattice basis and Fourier transform the spin- and bo- 
son operators on each sublattice separately: for sites r^ 
belonging to the A-sublattice we define 






k 



k 



(2.29) 
(2.30) 



and for sites r%- belonging to the B-sublattice, 



s t ] = JjrT,'""' 1 ®*' ( 2 - 31 ) 



h = 



£e ifer >B fe , 



(2.32) 



where the wave-vector sums are over the reduced (anti- 
ferromagnetic) Brillouin zone. The quadratic part H 2 = 
H^ + H^ of our effective boson Hamiltonian becomes 

H 2 = J n Sy2[Z h (A{A k + BiB k ) 



k 



tnt 



+m 2 7k (BlA k + A{B k )] , (2.33) 
where 7*, = J k /Jo with 



Note that 



N 



J Q S = 2DJS = ho/2. 



(2.34) 



(2.35) 



To completely diagonalize H2 we first introduce the sym- 
metric and antisymmetric combinations 



Ct 



1 
72 



[A k + aB k ] 



a = ±l , 



(2.36) 



and then perform a Bogoliubov transformation, 



with 






kn- 



Uka -CTVka 
-aVka u ka 



* 
* 



ka 



— ka 



where 



Ukcr = 



Vka 



Z h + am 2 ~f k + e ka 



2eka 



Z h + am 2 -f k - e ko 



2eka 



(2.37) 

(2.38a) 
(2.38b) 



with 



£fc(T = 



(Z h + crm 2 ~/k) 2 - («- 2 7fe) 2 
[Z h + ajk] [Zh - cr«7fe] 



1/2 
1/2 



(2.39) 



that 










U lo+ V ko = 


Z h + am 2 j k 

tko 


(2.40) 




^UkaVka - 


n 2 j k 


(2.41) 



Efeo 

Within linear spin- wave theory 5h — and hence Zh = 1 , 
but the factor Z^ will deviate from unity if we take higher 
orders in 1/ S into account. Since the above transforma- 
tions are canonical, our magnon operators ^ko satisfy 
the usual bosonic commutation relations, 



[*fc<r. **,'*'] =Sfe,fe"W . 



(2.42) 



In terms of the new operators ^ko the quadratic spin- 
wave Hamiltonian H2 is diagonal, 



H2 — 2_^ Eko 



>t 



1 



k <y 



*L*>™ + 5 



E, 



(i) 



with the magnon dispersions 

Eko — JoSeko 
The constant 



E, 



(i) 



N 



-,2^h 



2 -Z h J S = -NDJS Z s 



(2.43) 



(2.44) 



(2.45) 



is the l/S-correction to the ground state energy due to 
longitudinal spin fluctuations. The total l/S'-correction 
to the ground state energy is obtained by adding the 
zero-point energy of the transverse spin- waves to E^ 



(i) 



E, 



(i) 



E, 



(i) 
oil 



oE £fe » 



-NDJS 2 



Ci(ft) 
S '' 



(2.46) 



1 






Ci(h) = -^(Z h -e k o). 

ha 



(2.47) 



In the long-wavelength limit we obtain to linear order in 
5h = h — h c m and to quadratic order in fc, 



fe+ = mh c h + C.fc , 
2 _ „2„> xu , „2 .2 



(2.48a) 
E'i_ -- n z mh c 8h + ctk z . (2.48b) 

For small m the spin-wave velocities are 

c 2 + = C 2(l-3m 2 ), (2.49a) 

c 2 _ = cl(n 2 + 2m 3 5h/h c ), (2.49b) 



where cq is the leading large- S result for spin- wave ve- 
locity for h = 0, 



c = 2y/DJSa. 



(2.50) 



At the level of linear spin-wave theory we may approxi- 
mate the canting angle by its classical value $o> which is 
determined by the condition Sh — 0, or equivalently 



m = sini^o = h/h c . 



(2.51) 



This result can also be obtained by minimizing the clas- 
sical energy Eq in Eq. (|2.4[) . The gap of the dispersion 
Ek+ is then simply given by h, while the dispersion Ek- 
is gapless with spin-wave velocity 



c_ = c n = Cm/1 - 



h 2 
h 2 ' 



(2.52) 



C. Hermitian field operators 

In the usual l/.S'-approach one now substitutes the re- 
lations between the original Holstcin-Primakoff bosons 
bi and the magnon-operators ^ko into Eqs. (|2.22[ 12.241 
12.271 I2.28[) . This yields rather lengthy expressions in- 
volving momentum dependent vertices. However, if one 
is only interested in the transverse staggered spin fluctua- 
tions, it is better perform another transformation which 
separates the staggered from the uniform spin fluctua- 
tions. Therefore we express the magnon operators ^ko 
in terms of two Hermitian field operators Xko and Pko 
achieving the natural normalization on a lattice as fol- 



^ko = Po 



Vkc, 



Xko 



\J2vka 



--Pk 



(2.53) 



where the phase factors p + = —i and p_ = 1 are chosen 
for later convenience. Here the dimensionless factors i/ko 
are defined by 



Vko 



Eko 



±ko 



(2.54) 



where 



D. Spin- wave interactions 



A fecr = 2J Sz kcr = h c z ka , (2.55) 

and 

Zka = [Uka + V ka ] 2 e ka /2 

= [Z h + {n 2 + <jm 2 ) lk ]/2. (2.56) 

Note that Zh = 1 to leading order in 1/5, so that to this 
order 



where a = n 2 — 



Zk+ 

z k - 

9 

in 



(l + 7fc)/2, 
(l + 07fc)/2, 



(2.57a) 
(2.57b) 



In particular, for k — ► we have 



z k+ — ► 1 and Zfc_ — > (1 + a)/2 = n 2 . One easily verifies 
the canonical commutation relations, 



[X ka ,P k >cr>} = id, 



k.—k'Oa.a 



(2.58) 



The quadratic part of the spin-wave Hamiltonian can 
then be written as 



Ho 



1 



E A ^ 



ka 



P-koPko + V ka X- k<J X ko 



E 



(i) 



o|| • 
(2.59) 



In contrast to the lattice normalization of Eq. (|2.53[) in 
Ref. [l6| we focused on the continuum limit to exhibit 
the relation with the NLSM. In that case a continuum 
normalization of the fields is more convenient, 



** 



Pa 



Xo 



2VE k0 



E k<J H-k 



iXv l ®ka 



(2.60) 



where xo — {2Joa D )~ 1 is the large-5 limit of the uniform 
transverse susceptibility for h — 0. The continuum fields 
fulfill the commutation relation 



\Hka,®k'a>] = iV8h,-k'5tr,<T' 



(2.61) 



The relation between lattice and continuum normaliza- 
tions is 



n fco 



4>i 



,D 



N 
Sz kc 



-X kai 



\/NSz k<y P kt7 . 



(2.62) 
(2.63) 



Our spin-wave Hamiltonian (|2.43|) in continuum normal- 
ization can be written as 



H' 2 = ^zZ^o^-^ka+XaEl^^tlk,, 

ka 



t" ^o|| ■ 
(2.64) 



The field Ii ka corresponds precisely to the continuum 
field representing transverse staggered spin fluctuations 
in the non-linear sigma model 9 -. However, here we would 
like to calculate also short-wavelength properties on a lat- 
tice, so that we shall work with the lattice normalization 
(12351) . 



In order carry out the 1/5-expansion using the opera- 
tors X ka and P kcr defined in Eq. (|2.53p . we should first 
express the interaction part of the bosonized Hamilto- 
nian in terms of these operators. To obtain the leading 
1/5-correction to linear spin-wave theory, it is sufficient 
to approximate the effective bosonized Hamiltonian by 



H^E c Cl l + H 1 +H 2 + H 3 + H 4 



HI 



(2.65) 



where H± = H\ + H^. Later we shall use the phase 
space path integral to derive the effective action for 
staggered fluctuations. All expressions in the Hamilto- 
nian should therefore be symmetrized whenever powers 
of non-commutating operators are encountere d 18 ! 19 ' 20 . 
Only after symmetrization we may replace the field op- 
erators by numbers. If A\A 2 ■ ■ ■ A n is a product of op- 
erators consisting of X ka or P ka in arbitrary order, the 
symmetrized product is 



{iii 2 ---i„} = - T ^ip 1 ip 2 



■A, 



(2.66) 



where the sum is over all n\ permutations of 1, . . . , n. We 
obtain from Eq. (J2.27P for the linear part of the Hamil- 
tonian, 



H x = nShVSNPo-. 
The part H 3 in Eq. (]2.28[) can be written as 



(2.67) 



fffl 



N hrX 



85 



Po- 



22 ^fcl+fc: 



2+k 3 ,0 



k 1 k 2 k 3 



-r pxx (fc i; fc 2 ,fc 3 ){p fcl -i%_x fc3 _} 



1 



-y T -+i (ki;k 2 ,k 3 )P kl _x k2+ x k3+ 

-T pxx (k i; k 2 ; k 3 ){P kl+ X k2+ }X ks _ 
~^-++(kf,k 2 , k 3 )P kl _P k2+ P k3+ 



¥ 



ppp 



(fcl, fe 2 , fe 3 )P fcl _P fc2 _P fc 3 



(2.68) 



where the vertices are 

r pxx (k i; k 2 ,k 3 ) = 



h r X 



^PXX 



(ki;k 2 ,k 3 ) 



7 85 



T p f x (k i; k 2 ;k 3 ) 

r p ^(fe i; fe 2 ,fc 3 ) 

rf^ffci.fca.fca) 



'8S 

hr.X 



V8S 

h r X 



'8S 

h r X 



7fei. 


(2.6ya) 


[7fei lk 2 7fc 3 ] . 


(2.69b) 


lk 2 , 


(2.69c) 


7fei) 


(2.69d) 


[7fei + 7fe 2 + 7fc 3 ] • 


(2.69e) 



Explicitly, the symmetrized products in Eq. (I2.68P are 



{P1X2} 



Pi,Xi 



(2.70) 



{-Pl^2-^3/ — o (P1X2X3 + X2X3P1) 

+ i(x 2 A^3 + x 3 A^2) 

6 



Pi, X2X3 



(2.71) 



where [^1,^2]+ = A1A2 + A2A1 is the anti-commutator 
and we have abbreviated P kl by Pi and analogously for 
the other labels. 

Finally, consider the part H4 — HJ + H^ of the Hamil- 
tonian involving four boson operators, which according 
to Eqs. (j2~2"2")l and (j!T2"4"l) is given by 



H 4 



y^ Jij l TiiUj + - [mbibj + biUjbj 

ij L 

+b\b) nj + b\n t b]] 1 



'\ Yl J V \ Uin 3 ~ 1 [ n ^ h ] 



b,tf 



j n J 



+b\riibj + b\rijbj 



(2.72) 



Expressing H± in terms of the operators P ka and X ka 
defined in Eq. I|2.53[) and symmetrizing all expressions 
containing non-commuting operators we obtain 



Ha = E, 



(2) 

Oil 



SH 2 + H' 4 , 



where 



E, 



(2) 



NDJS 2 a 
(25) 2 



(2.73) 



(2.74) 



is a l/S^-correction to the classical ground state energy, 
and 



5H' 2 = \Y. \r*{k)P- ka P ka + TX{k)X_ ka X k!y 

fcfT 



is a l/S'-correction to H2* The vertices are 
rf(fc) = ^-a{l + a lk ), 



r?(fe) 



AS 
h c 

45 



(a - cr7 fc ). 



(2.75) 

(2.76a) 
(2.76b) 



Finally, the properly symmetrized quartic part H' A of our 
spin-wave Hamiltonian is given in Appendix A. For our 
purpose it is only important that the corresponding inter- 
action vertices are non-singular functions of the external 
momenta and are analytic functions of h 2 . 



III. EFFECTIVE ACTION FOR THE 
STAGGERED SPIN FLUCTUATIONS 

In Ref. [7j the precise relation between the magnon 
quasi-particle operators of the l/S'-expansion and the 
continuum fields H ka representing transverse fluctuations 
of the staggered magnetization has been established. In 
this section we shall use this relation to derive the ef- 
fective action for the staggered spin fluctuations for the 
Hamiltonian (|1.1|) retaining sub- leading l/S-corrections 
and short wave length fluctuations in the entire Brillouin 
zone. 

For weak magnetic fields, the operators P a corre- 
spond to transverse fluctuations of the total spin, while 
X a describe staggered (antiferromagnetic) spin fluctua- 
tions. To calculate the self-energy of antiferromagnetic 
magnons, we can therefore eliminate the degrees of free- 
dom associated with the generalized momenta P a . This 
is most conveniently done using path integration. The 
appropriate path integral in our case is the imaginary 
time phase space path integra l 18 ' 19 . Recall that for a 
one-dimensional quantum mechanical system with posi- 
tion operator X, momentum operator P, and Hamilto- 
nian H(P,X) the partition function can be written as 



2= J V[P,X]cxp I I dr 



r)X 

iP^-H s (P,X) 



(3.1) 
where H S (P,X) is obtained from the Hamiltonian 
H(P,X) by first symmetrizing H(P,X) with respect 
to the ordering of the operators X and P, and then 
replacing the operators by their eigenvalues. In prin- 
ciple, ambiguities associated with the operator order- 
ing in the phase space path integral can always be re- 
solved by going back to the discretized definition of the 
path integra l 18 ' 19 . However, recently Gollisch and Wet- 
terich 20 ' 25 showed that in the continuum notation the 
symmetrization prescription leads to the same result as 
the more fundamental discretized definition of the phase 
space path integral. The Euclidean action corresponding 
to our spin- wave Hamiltonian is of the form 



S[P a ,X a ] = YSi[P a ,X, 



2=0 



(3.2) 



where Si\P a ,X a \ contains I powers of the fields. To ob- 
tain the effective action S e s[Xa] for the staggered fluc- 
tuations, we integrate over the generalized momenta, 



s*[x r ] = r V [p a ] e s[p„ 



.x„ 



(3.3) 



Within the Gaussian approximation (corresponding to 
linear spin- wave theory) we truncate the expansion (|3.2[) 
at the term I = 2. The relevant contributions to 
S[P a ,X a ] can be written as 

So = 0[E* + E§;h 



Si[P-] = PnShy 



(3.4) 
(3.5) 
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and 



S 2 [Pa,X a ] = ^YX^a{P-KaPKa+vLX^ Ka X Ka ) 



K,a 



-W{P-Ko-X Ka — X_ KrJ P Ka ) , (3.6) 



where the last term in Eq. (|3.6|) corresponds to the mea- 
sure term iPdX/dr in the phase space functional integral 
(|3.1|) . The fields Pj£a and Xko are defined by replacing 
the operators Pka and Xka by quantum fields Pka(j) 
and X^o- ( T ) depending on imaginary time r and expand- 
ing the fields in frequency space, 



Pka(T) = J2 e ~ iUTP > 



Ka, 



Ka- 



(3.7a) 
(3.7b) 



We combine momenta k and bosonic Matsubara frequen- 
cies iu to form a composite label K = (fc, ioj). In general 
the canting angle can be determined from the condition 
that the functional average of the field Pk=o,- vanishes, 



(Po-) = 0. 



(3.8) 



Eq. (|3.8|) defines the correction 5h = h — h c m = h — 
h c sin ■& and hence the sine of the renormalizcd canting 
angle sin # = m = (h— Sh)/h c . Within the Gaussian ap- 
proximation this implies 8h = 0, leading to the classical 
result (|2.5ip . Hence 5i[P_] = within this approxima- 
tion and the effective action for the fields X a is given by 
the Gaussian integral 



-Smtt[X a 



,-S 



V[P a ]e 



— £>2 [P<r iX-a 



(3.9) 



Carrying out the integration, we obtain in Gaussian ap- 
proximation S c s [Xa] = So + 5g ff [X a ] , where 



-vi -.i = ?E sL+ " 2 



M r 



2 4r^ Afco- 

Ka 



X_ Ka X K<7 . (3.10) 



At long wavelengths this action has the same form as the 
corresponding Gaussian part of the action of the NLSM. 
However, in contrast to the NLSM, our action is defined 
on the lattice so that fluctuations on all wavelengths are 
included. The Gaussian propagator of the X CT -field is 
thus 



(Xko-Xk'o-')o — dK,-K>8o-o-'{0Ako-) 



EL+^ 



(3.11) 



The other propagators are within Gaussian approxima- 
tion 

E 2 
(PkvPk'<t>)o = 6 K ,-K'S^(f3A ka y 1 k ° , ,(3.12) 

^ka + UJ 

(X K oPk>o>)o = Sie-K'S^tfAto)- 1 „ 2 kaU 2 .(3.13) 



Here the symbol (. . .)o denotes functional averaging with 
the Gaussian action 5*2 \P a , X a ] . Note that the formal 
sum J2u>(XKo-P-Ka)o represents the expectation value 
of the symmetric operator ({XkaPka})o = 0, so that 
we should regularize formally divergent Matsubara sums 
using a symmetric convergence factor cos(cj0 + ), 



({XkaPka})l 



3(W0+) 



(3^ Ef 



0. 



(3.14) 



ka 



The higher l/S-corrections to S e ff[X a ], can now be ob- 
tained by including the spin-wave interactions perturba- 
tively. Therefore we rewrite Eq. (|3.3[) as 



(3.15) 



<Scff [A CT ] — Sq + S cS [Xa] + Srff [Xa] 



where the interaction part S£g [X a ] is defined via the fol- 
lowing functional average, 



S%[X„] = -lnw 
= -In 



JV[P]e- s ^ p ^ x ^e 



P 

S\ n k [Pa-)X a 



fV[P]e- s ^ p °< x ° 



(3.16) 



where 



S int [P a ,X a ] = SilP-l + ^SilP^X,]. (3.17) 

The leading correction of relative order 1/S arises from 
the first order correction due to Si[P a , X a ] corresponding 
to H 4 defined in Eqs. (j2~73l [27741 [27751 [Pj) . and the 
second order corrections due to the sum of Si[P_] and 
S 3 [P a , Xa], corresponding to H' « H x +H z in Eqs. (J27B7| 
and (|2.68[) . Note that to order 1/S the difference Sh = 
h — h c m and hence S*i[P_] are finite, so that the condition 
for the renormalized canting angle reduces to 



Pb-(Si[Pb-] + Sa[P„X ])) =0 



(3.18) 



Performing the Gaussian averages we obtain to first order 
in 1/5, 



8h = m[l-C 2 (h)]^, 



(3.19) 



with the numerical constant 

C%{h) = -TT^iula+V^-ajkiuka+CrVka) 2 ] 



kfj 



1 v-> 1 - 7fe - cr^ 2 7fc 



N £—1 



N 



Efeo 



(3.20) 



Our condition (|3.19[) leads to the same l/S'-corrections 
for the canting angle as in Ref. [2l[ and thus yields the 
same result for the uniform magnetization. Note that 
Si [P-] is of order 5 -1 / 2 and should be taken into account 



on the same footing with S 3 [P a , X a ] in second order per- 
turbation theory to collect all corrections of relative order 
1/5. Using Eq. (|3.19j) we obtain for the total contribution 
of order S* -1 / 2 to the action S^Pj, A CT ] corresponding to 
H' in Eq. (l2"TTD . 

S'[P a ,X a ] » SilP^ + SaiP^X.] 

FN krX 



-(3 



2 V8S 



C 2 (h)P - 



W£ E ^ 



+_ftT 2 +if 3 ,0 



K 1 K 2 K 3 



2^r (fei;fe 2 , k 3 )P Kl ^X K2 -X K3 ^ 

-^T p ^(k i; k 2 ,k 3 )P Kl ^X K2+ X K3+ 
-r™(fe i; fe 2 ; fc 3 )P Kl+ A X2+ A K3 _ 

~r^(fe 1 ;fe 2 ,fe 3 )Px 1 -Px 2 +Pif3+ 

~rf^(fci,fc2,fc3)P^_Pif a -Pif,- .(8.21 ) 



The leading correction to the Gaussian approximation 
for the effective action ScslX,,] is of order l/\/S, 



7(1/2) 



SiTTO = (S'[P^X a ]) P , 



(3.22) 



where the subscript indicates the power of 1/5. The 
1/5-correction is 



S^[X a ] = (S 4 [P a ,X a ]) P 
^((S'[P a ,X a }-(S'[P a ,X a }) P f /p 



(3.23) 



To calculate the Gaussian average in Eq. (|3.22j) we use 
the fact that averaging the field Pko for fixed X yields 



(P, 



Kcr/P 



A 



-X 



Ko 



k n 



(3.24) 



After proper symmetrization of the vertices we obtain 



o(l/2)ry 



"V at E ** 



+^2+^3,0 



K!K 2 K 3 



' -f { ll-{K x ,K 2 , K 3 )X Kl _X K2 _X K3 

-^l + (K i; K 2 ,K 3 )X Kl -X K2+ X K3+ \ (8.LT.) 



with 



i(3) 



rL°i_(^i,^ 2 ,ir 3 ) 



/i^ A 



Tfci^i 7fc 2 ^2 7fe3 w 3 



V8S"LA fcl _ Afc 2 - ' Afc 
(7fei + 7fc 2 +7fc 2 ) w i w 2W 3 



A fel _A fc2 _A* 



, (3-26) 



v_ 



A-1 



V- 



+ 




#2 



#1 



- -. ^3 



FIG. 2: (Color online) Graphical representation of the inter- 
action vertices V+(Ki, K2, K3) and V-(K\,Ki,K-i) defined 
in Eqs. (I3.28|l and (I3.29|l . Solid lines represent the gap less 
field X-, while dashed lines correspond to the gapped field 
X+. The shape of the symbols reflects the symmetry of the 
vertices with respect to the permutation of the labels. 



T^ ++ (K i; K 2 ,K 3 ) 



h r X 



7fe 2 ^3 . 7fc 3 ^2 . 

~~ 1 1 — 1 r 



(7fei - 7fc 2 " 

7/01^1^2^3 



7fc 3 ) 



*fc 3 + 



*fc 2 + 



*fc 2 + z 



A fcl _ 

(3.27) 



Actually, the terms cubic in the frequencies which are 
due to the cubic terms in the Pk<j in Eq. (|3.2ip can be 
omitted, because the contribution of these terms to the 
self-energy of the A-fields is frequency-independent to 
order 1/S. Since we are only interested in the frequency 
dependent part of the self-energy, we may thus replace 



r^i_{K u K 2 ,K a 



V-(K l3 K 2 ,K 3 ) 



h r X 



7fci w i . 7fe 2 ^2 . 7fc 3 ^3 



^(3) 



8S*L^ fcl 

_(K i; K 2 ,K 3 ) 



hc.X 



ilk, -7fe 2 ~7fc 3 ) 



^fc 3 - 
V+(K U K 2 ,K 3 ) 

u\ 7fc 2 ^3 



(3.28) 



7fc 3 ^2 



A, 



Mt3 + 



^2 + 



(3.29) 



Graphical representations of the interaction vertices 
V a {K 1 ,K 2 ,K 3 ) are shown in Fig. [5] 

At this point we can make contact with the NLSM, 
which is an effective low-energy theory for staggered spin 
fluctuations. In the presence of a uniform magnetic field 
the Euclidean action of the NLSM isiS^ 3 . 



Snlsm[^] = 



I f dr f d D r[£(d^r 

+c- 2 (d T n-ih x fl) 2 }, (3.30) 



where the unit vector Jl(r, r) represents the slowly fluctu- 
ating staggered magnetization, p s and c are the spin stiff- 
ness and the spin- wave velocity at temperature T = 0, 
and dfj = d/dr^ is the spatial derivative in direction 
fi = 1,...,D. The model (|3.30[) can be obtained from 
the corresponding NLSM for h = by substituting 
d T — ► d T — ihx. Although this procedure does not explic- 
itly take into account the magnetic field dependence of 
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the spin-wave velocity and the spin stiffness, one usually 
argues that c and p s in Eq. (|3.30[) are effective parame- 
ters, implicitly including the effect of the magnetic field. 
However, this procedure is based on the assumption that 
in the presence of a magnetic field the magnon dispersions 
can be characterized by a single spin- wave velocity c(h). 
From Eqs. (|2.49ap and (|2.49b|) it is clear that this assump- 
tion is not justified, because the dispersion of spin-wave 
mode polarized parallel to the magnetic field involves a 
different spin- wave velocity than the mode polarized per- 
pendicular to the magnetic field 16 . Apparently, there are 
no published calculations of the l/S'-corrections to the 
magnetic field dependence of the spin-wave velocity. In 
the following section we shall show that in dimensions 
D < 3 the magnetic field dependence of the spin- wave ve- 
locity c_ (h) of the gapless magnon mode is non-analytic 
in h 2 . 

To make contact with our spin-wave approach, let us 
consider the interaction vertex due to the magnetic field 
in the NLSM. Therefore we rewrite Eq. ([330]) as 



Aka- ~ h c , so that we obtain 



r^L (Ki,K 2 ,K 3 ) 



[u>i + uo 2 + UJ3] = 0, 



T% + (K i; K 2 ,K 3 ) 



[— Wl +LU2+ W3] 



-2-=wi, 

V8S 



(3.33) 



(3.34) 



where we have used the fact that loi + lo 2 + lu 3 = by 
energy conservation. Finally, using the relation (|2.62|) 
between continuum and lattice normalization of the field 
representing the staggered spin fluctuations, it is easy to 
see that for weak magnetic field the continuum limit of 

our lattice action S^ s [X a ] in Eq. (|3.25[) reduces to the 
cubic term in the expansion (|3.32p of the NLSM. 



IV. FREQUENCY DEPENDENT PART OF THE 
SELF-ENERGY TO ORDER 1/5 



'NLSM 



[fi] 







dr Id 



iD 



D 

r[^(^fi) 2 + c- 2 (9 T n) 2 



pv\ti> + \ 



I dr J d D r{h-to) 2 
I dr I d D rM ■ (to x d T to) , (3.31) 



where \ = Ps/c 2 and M = \h. Choosing the coordinate 
system such that the staggered magnetization points in 
direction e z and keeping in mind that h = he x , we now 
set to — \/l — n 2 e z + II and expand Eq. ()3.31j) in powers 
of the transverse fluctuations II. Retaining only terms up 
to cubic order in the fluctuations II = U + e x + H-e y we 
obtain in momentum-frequency space, 



-PV^-h 2 
2 



Snlsm [to] fH 

1 Jk a 
-i X h I dr f d D rU 2 + d T n^+0(U 4 a ), (3.32) 



where m 2 _ — and m 2 ^ = h 2 . At the first sight, the cu- 
bic interaction in Eq. (|3.32p does not resemble the cubic 

term S^ /2) [X a ] in Eqs. (|3.25H3.27|I . However, the NLSM 
is only valid to leading order in the derivatives, so that for 
a comparison with Eq. (|3.32p we should expand the ver- 
tices (J3.26I) and (|3.27p to leading order in momenta and 
frequencies. Moreover, for small h we may approximate 



Defining the non-interacting propagators of the stag- 
gered spin fluctuations, 



GoA K ) = 



A 



krr 



EL+" 2 



(4.1) 



and expressing the corresponding interacting propagators 
in terms of the self-energies Yl a (K), 



G- 1 (K) = G^ 1 a (K) + ^(K), 



(4.2) 



the leading frequency dependent contribution to the self- 
energy correction of the gapless magnon mode can be 
written as 

£_(*■) = ^tT.T.^Ak^GoAk' + k) 

P K' a 

xV 2 (K,K',-K -K'), (4.3) 

while the self-energy of the gapped magnon mode is 



K' 



xVl(K',K,-K-K') 



(4.4) 



where we have used V a (— K, — K', K + K 1 ) = 
— V a (K, K' , —K — K 1 ). The corresponding Feynman dia- 
grams are shown in Fig.[3l The frequency integrations in 
Eqs. (|4.3p and (|4.4p can now be performed analytically; 
the relevant integrals are 



I (n) (Bi.^.w) = 



E[ l - 



dx 



_ 00 2nlx 2 + E 2 }l(x + u;yi + E 2 } 
(E 2 + iuj) n 



E 2 2 



{Ei 



E 2 [E 2 - (E 2 + iu) 



, (4.5) 
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S_ = -i 



s + = -^ 







FIG. 3: (Color online) Feynman diagrams of the self-energy 
corrections to second order in the three-legged vertices, see 
Eqs. (|4,3[l and l|4.4|l . The slashed tadpole diagrams give 
frequency-independent contributions of order 1/S which are 
analytic functions of the magnetic field. Since in this work 
we are only interested in the frequency dependent part of the 
self-energy, we shall omit the tadpole diagrams. 



J 



where n = 0, 1, 2. Explicitly, 



j(o) 

j(D 

j(2) 



Ei + E 2 



2E 1 E 2 [{E 1 + E 2 ) 2 + uj 2 

LO 

~ 2E 2 [(E 1 + E 2 ) 2 + u 2 } 

E 2 {E X + E 2 ) + LO 2 
2E 2 [{E 1 + E 2 ) 2 + to 2 ]' 



(4.6a) 
(4.6b) 

(4.6c) 



The result for the self-energies can be written as 



E-W 



£+(#) 



h\\ 2 2 

16S N 



n 



Zq-\-Zk—q-\- 



Zq—Zk—q- 



M 2 (fe, g)li°i(tw, fc, q) + 2M (fc, q)M+(k, q)P+i(iu, fc, q) 



-Ml(k,q)I ( +l(iu;,k,q) 
r(°), 



r(i) 



fegA 2 2 

16S AT 



M 2 (fc, q)I^(ioj, k, q) + 2M_(fe, g)M_(g, k)I^(iuj, k, q) 
hMl(q,k)I { y_{iuj,k,q)\\ , 

Y, z q -z k - q + \Ml(q, k)I ( °l(iu, k, q) + 2M+{q, k)M (q, k)I ( y + (ico, fc, g) 



+M 2 (<7,fc)/i 2 j^,fc,q) 
I 



(4.7) 



(4.8) 



where 



I%{iu,k,q) =u> 2 - n I^(E qa ,E k _ qlT ,,oj), (4.9) 



and we have introduced the functions 



M (k,q) 
M+(k,q) 
M_(k,q) 



lq 



lk 1q Ik — q 



%k—q-\- %k- 

lq _ lk-q 

1 

%k—q-\- Zq-\- 

lk lk-q 

%k— ^k—q — 



, (4.10a) 
(4.10b) 
(4.10c) 



For later reference we note that 



M (0,q) = ^_ + 2 -^l, 

Zq+ Z - 

M+(0,q) = 0, 
M_(0,q) = — -^-, 

ZQ- Z q - 

M o (fc,0) = — + — , 

Zk+ Zk- 



Af + (fc,0) = 



1 



7fe 



z k+ Zq + 



M_(fc,0) = 0. 
Furthermore, if both fc and q are small 



(4.11a) 
(4.11b) 
(4.11c) 

(4.11d) 

(4.11e) 

(4.11f) 



M_(fc,q) 



4Dn 4 



[q 2 -2k-q] +0(k\q 4 ,k 2 q 2 



(4.12) 
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V. RENORMALIZATION OF THE GAPLESS 

MAGNON 

A. Spin-wave velocity 

We now show that in dimensions D < 3 the leading 
1/5-correction to the spin- wave velocity c—{h) of the gap- 
less magnon is non-analytic in h 2 . Therefore we expand 
for small w and \k\, 

S_(fe,zw) = f uj 2 + jiu) 2 , 

+ } 2 uj a + huo 2 uj 2 k + Uwi + 0(^ 6 ), (5.1) 

where u> and toy- = c_|fe| are assumed to have the same or- 
der of magnitude and c_ = cqu 2 is the spin-wave velocity 
within linear spin- wave theory, see Eq. (|2.52| . To calcu- 
late the renormalized spin- wave velocity we may neglect 
in Eq. (|5.1[) the terms of order cu 4 involving the coeffi- 
cient / 2 , /3 and /4. Using Eqs. (14. ip and (|4.2p we obtain 
for the infrared behavior of the propagator of the gapless 
mode 

Z-h r n 2 



G-{k,iu) 



u* 



c 2 k 2 ' 



(5.2) 



Introducing the dimensionless constants Fq and F\ , 

F = h c n 2 f , F x = h c n 2 h, (5.3) 

the wave-function renormalization factor Z_ can be writ- 
ten as 



Z- 



1 



1-Fn 



1 + i^o 
and the renormalized spin-wave velocity c_ obeys 

1 + Fi 



(I 



1 



1 + F-l-Fq. 



(5.4) 



(5.5) 



The constants fo and / 2 associated with the expansion in 
powers of frequencies for vanishing external momentum 
can be obtained by expanding S_(fc = 0, to) in powers of 
lu 2 . Using Eq. $77} and Eqs. (|4.11aft4~TTf]) one gets 



E_(0,m/) 



h 2 X 2 2 v^f o ,,2, n x r (o) 



I^^2.(^ + M o(0,q)/n(-,0, 9 ) 

z2_M 2 (0,q)/i°l(^,0,q)} 



2 h 2 c X 2 2 



[7^ + ^(27,-1) 



165 JV^ J. £ q+ [(2£ q+ ) 2 +^ 2 
fc-7,] 2 



£ <J _[(2 J B <2 _)2+ W 2] f '"■" , 

Using h 2 X 2 — 4n 2 h 2 , we obtain for the first two coeffi- 
cients in the frequency expansion, 

n 2 h 2 2^([ lq + 1 -^(2 lq -lf 



k 



165 N 



E 



[fa 



^ 



7 9 J 



E\- 



(5.7) 



/a = - 



n 2 /i 2 2 v^ 



165 N ^ 



w 



^p-^ Z 7<3 



1)] : 



4B 5 



4^_ 



(5.8) 



Keeping in mind that z q _/zo- — 7q = 0(q 2 ) for small 
q, it is easy to see that in the domain of small magnetic 
field (h <C /i c ) the integrals on the right-hand sides of the 
equations above are dominated by the first term involving 
the gapped mode E q+ . More precisely, the relevant ul- 
traviolet cutoff for the momentum integrals in Eqs. (|5.7[ 
15. 8p is the inverse of the length scale 



£ = Co/ft. 



(5.9) 



In D < 3 the contribution from wave-vectors in the 
regime |q|£ < 1 gives rise to contributions to the magnon 
self-energy which are non-analytic in h 2 . Keeping in 
mind that for small field the magnetic length £ is large 
compared with the lattice spacing, we may calculate the 
leading non-analytic magnetic-field dependent contribu- 
tions to Eqs. (|5.7[ 15. 8p by expanding the integrand in 
powers of q. 

We find that the leading magnetic field dependence 
of the spin-wave velocity c_ associated with the gapless 
mode is determined by /q. Since we are only interested 
in the non-analytic /^-dependence, we may set u»l, In 
the thermodynamic limit we then obtain for the domi- 
nant contribution to Eq. (|5.7p . 



2„D 



k 



h 2 a 



d D q 1 
25 J (2^W, 



(5.10) 



Consistently neglecting terms which are analytic in ft 2 , 
we may ignore the magnetic field dependence of the non- 
interacting spin-wave velocities, c± fa cq — 2yDJSa, 
so that energy dispersions are approximated by E q - s=a 
co\q\ and E q+ « \Jh 2 + c^q 2 . Using h c = 2y/Dco/a we 
obtain from Eq. (|5.7p for the corresponding dimensionless 
coefficient for 1 < D < 3, 



,D-1 



^0 = h c n /o = a D - 



(5.11) 



where m = h/h c = ha/(2\/Dco) is the relevant dimen- 
sionless magnetic field [see Eq. (|2.51[) ]. and 



an =2 V ' K 



D 



1/m 



,-D-l 



dy 



o y [l + 2/ 2 ] 3 / 2 ' 



Here 



>1-D 



K D = 



ir D / 2 T(D/2) 



(5.12) 



(5.13) 



is the surface area of the Z?-dimensional unit sphere di- 
vided by (2ir) D . In D < 3 we may take the limit 
1/m — ► oo in ap, so that 



'D 

a D = I — 



D/2 r(3f>) 



(5.14) 
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In particular, a 2 = 2/ir. In D — 3 the integral 0:3 de- 
pends for small m logarithmically on the upper limit, 



a 3 ~a 3 ln(l/m), a 3 



6^ 



(5.15) 



It turns out that the coefficient F\ in front of the k 2 - 
correction to the self-energy is for small h proportional 
to h 2 , so that for h <C A the dominant magnetic-field 
dependence of the spin-wave velocity is due to the term 
Fq in Eq. (|5.5p . We thus obtain for the leading magnetic 
field dependence of the spin-wave velocity of the gapless 
magnon 



1-Fo 



1 



J_\h\ 
irS h r 



D 



TT 2 S h 2 Mlhl 



(5.16a) 
L> = 3, (5.16b) 



where we have neglected magnetic field independent 1/S- 
corrections. Recall that within linear spin-wave theory 
the velocity c_ of the gapless magnon is analytic in 
h 2 = h 2 m 2 \ from Eq. (|2.52[) we obtain c_ ra Co[l — m 2 /2] 
for small m. We conclude that in dimensions D < 3 the 
dominant magnetic field dependence of the spin-wave ve- 
locity of the gapless magnon is due to spin-wave inter- 
actions. In Appendix B we show that the non-analytic 
dependence on h 2 predicted by Eq. (|5.16aP can be recov- 
ered numerically from in the expression for the magnon 
self-energy given by Zhitomirsky and Chernyshev 10 . 



B. Quasiparticle residue 

In view of the fact that the magnetic field dependence 
of the spin-wave velocity of the gapless magnon is domi- 
nated by spin- wave interactions, it is reasonable to expect 
that also the higher coefficients in the expansion of the 
self-energy of the gapless magnon for small wave- vectors 
and frequencies exhibit some non-analytic dependence 
on the magnetic field. Consider first the renormalized 
magnon energies E k o, which can be defined by 



EL = E 2 ka + A fcCT ReS ff (fc, E k(y + tO). 
The expansion for small wave-vectors is 



El- = ctk 



1 + i_(fc)fc 2 + 0(k 4 



(5.17) 



(5.18) 



It is well knownS 2 - that only if the coefficient A- is pos- 
itive a gapless magnon with momentum k can sponta- 
neously decay into two magnons with momenta q and 
k — q. Within linear spin-wave theory we obtain from 
Eqs. (f2739|) and (|2~44| in D dimensions 



EL 



k u 



l + ^_(fc)fc 2 + (9(fc 4 ) 



(5.19) 



E, 



k+ 



= h 2 



c 2 + k 2 



1 + A+{k)k 2 + 0(k 4 ) , (5.20) 



with 



A-{k) 
A + (k) 



a 
"T 

a 2 

' 4 



1 -2m 2 



1 



L>(l-m 2 ) 3 

1 -2m 2 
£>(1 - 3m 2 ) 



-EK 


, (5.21) 


M J 




1" 


. (5.22) 



Obviously, for m <C 1 the coefficient A_(k) is negative 
for all directions fc, so that to this order in spin- wave the- 
ory the gapless magnon cannot spontaneously decay at 
long wave-lengths. For larger m the coefficient A-(k) 
decreases and eventually vanishes at a critical m sr {k) 
which depends on the direction k. From Eq. (|5.21[) 
it is easy to show that the direction where m*(k) as- 
sumes the smallest possible value is given by the diag- 
onal k x = ... = ko, and that the associated minimum 
is m* = h*/h c = 2/V7 « 0.76. For the special case 
D = 2 this result has been obtained previously by Zhit- 
omirsky and Chernyshe\*iJi, who examined the leading 
1/5-correction in the regime /i» < h < h c numerically. 

Apparently, the leading l/S'-correction in the limit 
of small magnetic fields m — h/h c -C 1 has not been 
explicitly analyzed in Ref. [l0| . In terms of the ex- 
pansion coefficients introduced in Eq. (|5.1[) we obtain 
A-(k) = A-(k) + 6A-, where the 1/5-correction is 



6A- = c 2 a h c [f 2 - h + h] 



(5.23) 



Let us consider first the contribution from the coefficient 
fi related to the ui -term in the expansion of the self- 
energy S_ (0, iu) for small frequencies. Because for small 
h the integral defining f 2 in Eq. (|5.8[) the dominated by 
wave- vectors \q\ < h/co, we may approximate 



1„D 



h 



h 2 a 



85' 



d u q 1 

WWW 



(5.24) 



q+ 



The integral is easily evaluated to leading order for small 
m<l. Introducing the dimcnsionlcss coefficient 



Fo 



clh c f 2 



we obtain for D < 3, 

0D 



Fo 



D-3 



S 



+ 0(m D - 1 )], 



(5.25) 



(5.26) 



with the numerical coefficient 

{2VD^ D - 2 



ll 



D 



-K T 



(2VD 



D-2 



X D-1 

„ dx [i + x 2 Y>/ 2 
2 



30F 



r(- 



^»r«f: 



In particular, in two dimensions j3 2 — l/(487r). 
ously, for D < 3 the coefficient F 2 diverges for m 



(5.27) 

Obvi- 

■+ 0, so 
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that the contribution from the term f 2 to 8A_ is for suf- 
ficiently small m much larger than the linear spin-wave 
result (|5.21|) . It turns out, however, that the singular 
contribution to <L4_ due to fa is exactly canceled by a 
similar contribution from the coefficient / 3 . In order to 
extract the dominated contribution to / 3 , it is sufficient 
to approximate the magnon self-energy (|4.7[) by 



£_(fc,m>) 



2h 2 a D 



d D 



9 r(0) 



(27T) 



I^i(iu;, k,q). (5.28) 



D A ++ 



Expressing m — h/h c — ha/(2y'Dco) = a/(2\AD£) in 
terms of the length scale £ = cq/Ii associated with the 
magnetic field we may alternatively write 



l-Fn- 



/W£ 



3-D 






1 /a\D-i 



k 2 a 2 



&D - Pok 2 C 



(5.34) 



Expanding the right-hand side to second order in k and 
comparing with Eq. (|5.1|) . we obtain 



h 



h 2 a D 


f d D q 1 


3- 


lOcfo 2 ! 


85 J 



.(5.29) 



The integral can easily be carried out analytically with 
the result / 3 = / 2 + 0(m D_1 ). From Eq. (|4~7)1 we can 
also show that the term j\ is of order a 2 m D ~ 1 /S and 
can be neglected as compared with f 2 and / 3 . Because 
S A- involves the combination f 2 — / 3 , we conclude that 
the singular contributions proportional to m D ~ 3 cancel 
in 5 A- , so that the leading magnetic field dependence of 
A- is proportional to m D ~ l oc l/il" 0-1 . This is small com- 
pared with the linear spin-wave result but non-analytic 
in h 2 , similar to the leading magnetic field-dependence 
of the spin- wave velocity in Eqs. (|5.16a,l [5.16b|) . 

On the other hand, the singular magnetic field depen- 
dence appearing in the coefficients f 2 and / 3 does not 
cancel in the self-energy £_ (k, w + iO) off resonance. Re- 
taining only the singular contributions to Eq. (I4.7P we 
obtain with f 2 « / 3 



£_ (fe, u + iO) « -foiu 2 + f 2 Lo 2 (Lu 2 



c 2 _k 2 ) 



(5.30) 



where a D = a D {2^) 1 - D and D = f3 D {2^/Df- D . In 
particular, in D = 2 the leading momentum dependence 
of Zk- is proportional to fc 2 £a = k 2 coa/h. The higher 
powers in k become important for |fc£| > 1, so that the 
expansion (|5.34[) is limited to the regime \k\ < t;^ 1 <C a -1 
where the l/jS-correction is small compared with unity. 



C. Magnon damping 



Given the magnon self-energies T,,y(K) in Eqs. (|4.7I4.8|) 
and the renormalized magnon dispersions E k a, the 
magnon damping can be obtained from 



T ka = ^lu^ a {k,E ka +iO). (5.35) 

2-Efctr 



The corresponding renormalized magnon Green function 
for small m can be written as 



G-(k,iu) = Z_(iu>) 



h r n 2 



uj 2 + ~c 2 _k 2 ' 



(5.31) 



where the renormalized spin-wave velocity is given in 
Eqs. (|5.16al5.16b|) , and 



Z-(iui) 



1 



1 + F + h c n 2 f 2 L0 2 
1-F - (a 2 F 2 /c 2 )iu 2 . 



(5.32) 



After analytic continuation to real frequencies we obtain 
for the renormalized residue of the magnon peak for small 
m , 



= Z_(iuj^c-\k\) = 1-F 



F 2 k 2 a 2 



2„2 



1-F -(3 D 



k 2 a 



(5.33) 



Zhitomirsky and Chernyshev- have shown that in two 
dimensions one should self-consistently take into account 
the imaginary part of the magnon self-energy when eval- 
uating the integrals on the right-hand side of Eq. (|4.7[) . 
However, as long as we are not too close to the criti- 
cal field ft-*, the result for the magnon damping is non- 
singular even if we ignore the damping of intermediate 
magnons in Eq. (|4.7p . We therefore expect that a sim- 
plified version of Eq. (|5.35|) taking into account only the 
renormalization of the real part of the magnon dispersion 
yields a qualitatively correct estimate for the magnon 
damping away from /i*. 

To calculate the damping T k - of the gapless magnon 
for wave-vectors \k\ <C h/co = £ _1 , it is sufficient to 
retain in Eq. (|4.7[) only the terms involving the functions 

/__ (iu>, k,q), because the imaginary part of the functions 

I { +i(uj + iO, k, q) vanishes for w < 2h. Using Eq. (J4~T2|) 
we obtain for u> > 
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Im 



Ml(k, q)P_L{w + iO, k, q) + 2M_(fe, g)M_(g, k)lrL(w + z'O, fe, g) + M 2 (g, k)P_l{uj + iO, fe,g) 



(5.36) 



r 



where 






W(k,q) = 


\k-q\ q 


-2k- q) 2 




-2(k 2 -q 2 )(q 2 -2k-q). 


(5.37) 



Note that in the non-linear sigma model the contribution 
corresponding to Eq. (|5.36[) is neglected because the rel- 
evant vertex involving three gapless magnons is set equal 
to zero (see Eq. ([333])), which is correct to leading or- 
der in the derivatives. Hence, the damping of the gapless 
magnon cannot be obtained using the NLSM. To estimate 
the magnon damping we set lo = E k - and approximate 
the renormalized magnon dispersion by 



£ fe _«c_|fc|(l + I_fc 2 ), 



(5.38) 



where for simplicity we have replaced the direction- 
dependent coefficient A_(fe) defined in Eq. (|5.18|) by 
some angular average A- . At long wave-lengths we then 
obtain 



r fc - = 



D h 2 a D + 3 



d n q W{k,q) 



8(4£>) 2 S J (2ir) D 
x6(E k ^-E k - q --E q -). 



(5.39) 



As discussed in the textbook by Lifshitz and Pitaevskii 22 , 
in the long wave-length limit the energy conservation 
E k - = Ek-q- + E q - can only be satisfied for A_ > 0. 
From our discussion in Sec. IV Bl (see also Ref. [I(|) we 
know that this condition is only satisfied in a certain 
range h* < \h\ < h c of magnetic fields below the sat- 
uration field. We now restrict ourselves to this regime, 
without explicitly calculating the magnetic-field depen- 
dence of the coefficient A- > 0. If /i is not very close to 
the threshold fields ft* and ft c , we expect by dimensional 
analysis that A- /a 2 is a number of the order of unity. 
The energy conservation then implies that the allowed 
vectors q are almost parallel to the direction of k and 
satisfy q < k. In fact, it is easy to show that the angle 
f? between fc and q is $ ~ y / oAl(k — q) due to energy 
conservation so that for A_k 2 <C 1 we may approximate 



6{E k . 

and 



— Ek-q- — E„ 



5(0 - y/6Al(k - g)) 
^j6A-C-kq 



, (5.40) 



Ifc-qfl w k-q+— — (1-costf) 
k — q 

w {k-q){l + ZA_kq). 



(5.41) 



Keeping in mind that A-kq <C 1 we obtain from 
Eq. (IOTP . 



W(k,q) 



9kq(k-q). 



(5.42) 



The integrations in Eq. (|5.39[) are now elementary and 
we obtain for the damping of magnons with wave- vectors 
in the regime \k\ < H/cq ^C a -1 at zero temperature in 
D dimensions, 



Ek- 

where 



7£ ( h_ 
S \h c 



6A_ 



D-3 



,D+li. I2D-2 



(5.43) 



Id 



Kd-i / dx[x(l — x)] 
Jo 

-Kd-\2 



D-l 



6WD 



64V£> r(£>+i) 

In two dimensions we have 72 = 3/(128v / 2t'") and 



(5.44) 



r.- = |,- 



2 - 



\k 



3„3 



V6A- 



(5.45) 



The |fc| 3 -dependence of the magnon damping has been 
obtained previously by Zhitomirsky and Chcrnyshev^. 



VI. SUMMARY AND CONCLUSIONS 

The main result of this work is the discovery that in 
quantum Heisenberg antiferromagnets subject to a weak 
uniform external field the leading l/S'-correction to the 
self-energy of the gapless magnon is a non-analytic func- 
tion of h 2 in dimensions D < 3. We have explicitly 
calculated the leading magnetic field dependence of the 
spin-wave velocity and the momentum-dependent quasi- 
particle residue of the gapless magnon. At the first sight 
it is surprising that for quantum antiferromagnets in a 
uniform magnetic field at zero temperature the dimen- 
sion D — 3 plays the role of a critical dimension below 
which fluctuations lead to a non-analytic magnetic field 
dependence of the magnon spectrum. However, the gap- 
less magnons in our model can be viewed as an interact- 
ing Bose gas in the condensed phase 23 , where the Bogoli- 
ubov mean-field theory is know n 24 ' 25 to break down in 
dimensions D < 3. 
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Finally, let us point out that our hybrid approach be- 
tween l/iS-expansion and NLSM is a very convenient pa- 
rameterization of the spin-wave expansion, which should 
also be useful in other contexts. While the calculations 
presented here can (with some effort) also be carried 
out using the conventional parameterization of the 1/S- 
expansion, our hybrid approach greatly facilitates the 
identification of the frequency dependent contributions 
to the magnon self-energies which give rise to the domi- 
nant magnetic field dependent corrections to the magnon 



J 



spectrum. 
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APPENDIX A: QUARTIC SPIN- WAVE INTERACTION IN HERMITIAN FIELD PARAMETERIZATION 

In Hermitian field parameterization, the quartic part of the Hamiltonian H^ defined in Eqs. (J2.72I2.73P is 



H' 



N 



E ** 



+k 2 +k 3 +k 4 ,0 



felfc 2 fe 3 fc4 



— { fr++++(fci, Ai2, fe 3 , kA)Ph 1 +Ph 2 +Pk 3 +Ph i + +r^^.(ki,k2,k 3 ,k4)Ph 1 -Ph 2 -Ph 3 -Ph i 

+-jt (r (fci, &2, &3, k^x^-Xhz-Xhs-Xht- + r ++++ (fei, &2, fc3, k4)x kl+ x k2+ x k3+ x ki+ J 

+ ,„,,, 2 (r +H (fci, k 2 ; k-j,, k4)Pk 1 + Pk 2 + Pk 3 -Pk i - + T__ ++ (fci, k-z\ k-i, k4)Xk 1 -Xk 2 -Xk 3 + Xk i + 

+ (o\\2 \T++— (ki, k 2 ; k :i , k4)Pk 1 +Pk 2 +Xk 3 -X k4 - +T — ++ (fci, fc 2 ; k s , k4)P kl -P k2 -X k:i+ X ki+ 

+ 7^2 ( r ++++ (fei. fc 2; fe 3, fc 4 ){ Pk 1+ Pk 2+ Xk 3+ Xk i+ } + rffff(fci, fc 2 ; fcs, fc 4 ){ P kl - Pk 2 -Xk 3 -X kA -}) 



(2!) 2 - 
+r^5(fci;fe 2 ;fc 3 ;fc 4 ){P fel+ ^fe 2 + }{P fe3 -X fe4 _} 



(A.l) 



where the symmetrization symbol {. . .} is defined in Eq. (|2.66p and we have used 



{P1P2X3X4} = - Pi P2, ^3-^4 + -7 (^1+3, 0^2+4,0 + ^1+4, 0^2+3, o) 



(A.2) 



For convenience we now introduce the short notation 7^ = 71, jk 2 = 72 (and similarly for the other labels) and 
symmetrize the vertices whenever the interaction is symmetric with respect to the exchange of the field labels. For 
the vertices involving four fields of the same type we obtain 



F (fci, k2,ki, fci) 



r + +^ x (fci ,k 2 ,k 3 ,k4) 



r++++(fci,fc2;fc3,fc4) 



r^^(fei,fc 2 ,fc3,fe4) 



—^ 71 + 72 + 73 + 74 - 2q(7i +2 + 73+4) + (2 <-> 3) + (2 <-> 4 
166 L 

j^ [0(71 + 72 + 73 + 74) - 2«( 7 i +2 + 73+4) + (2 <-> 3) + (2 <-» 4)] , 
-^~ I -71 - 72 - 73 - 74 - 2q(7i +2 + 73+4) + (2 <-*• 3) + (2 <-> 4) I , 
-^ [-01(71 + 72 + 73 + 74) - 2a( 7 i +2 + 73+4) + (2 <-> 3) + (2 <-> 4) 



(A.3) 
(A.4) 
(A.5) 
(A.6) 
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The vertices involving two pairs of fields of the same type can be written as 



r55^ H X (fei,fe 2 ;fe3,fc4) = — ^|^7i +72 + 0(73 +74) - 2a(7i +2 +73+4)J, 

r++^(fci,fc 2 ;fc3,fc 4 ) = Yfiu 71 ~ 72 _a( ^ 73 + T4) -2^(71+2 + 73+4) 
r++55(fci,fc 2 ;fc3,fc 4 ) = Y^[ 3 (~7i -72 + 73 +74) 

— 2Q (7l + 2 + 73 + 4 — 7l + 3 — 72 + 4 — 72 + 3 — 7l+4) , 

r__++ (fci,fc 2 ;fc 3 ,fc 4 ) = Y7j^ [ 3a (~7i ^72 + 73 +74) 

—2a (71+2 + 73+4 — 71+3 — 72+4 — 72+3 — 71+4) , 
r+++f(fci,fc 2 ;fc 3 ,fc 4 ) = ^ [-7i -72 + a(73 + 74) -2a(7i +2 +73+4)], 
r (fei,fe 2 ;fe3,fc 4 ) = — -77; — 0(71 + 72) + 73 + 74 — 20(71+2 +73+4) 



And finally, there is one vertex without permutation symmetry connecting four different field types^, 
r++^(fei;fc 2 ;fe 3 ;fe4) = ^-|;[7i + Q (-72 + 73) -74 - 2a(7i +4 +72+3)J. 

I 



(A.7) 
(A.8) 

(A.9) 

(A.10) 
(A.ll) 
(A.12) 

(A.13) 



Note that the above vertices are analytic functions of the 
external momenta and of h 2 . On the other hand, if we 
express H^ in terms of the usual magnon creation and 
annihilation operators, we obtain vertices which are sin- 
gular for certain combinations of external moment a 5 ' 6 ' 16 . 



APPENDIX B: N UMER ICAL CONFIRMATION 
OF EQUATION (|5.l6a|) IN TWO DIMENSIONS 

In this appendix we briefly review the calculation of 
the l/S'-corrections to the field dependent spin- wave dis- 
persion in two dimensions as obtained within the conven- 
tional 1/5-expansion by Zhitomirsky and Chernyshev in 
Ref. [l(| • From the numerical analysis of this expression 
we quantitatively confirm our result given in Eq. (|5.16a[) 
for the linear magnetic field dependence of the spin- wave 
velocity associated with the gapless magnon. In our nota- 
tion the expression for the on-shell renormalized magnon 
energy E kcj given in Ref. [10| can be written as 



where the self-energy has the form 



-il/S 



(fe, iw) 



-,1/S, 



il/S, 



+zll s (k) + i: 1 4 l s (k). (b.2) 



The frequency dependent contributions to the self-energy 
are given by 



-a/s 



(fe, iw) 



Xll S (k,iw) = 



h 2 c \ 2 2 ^ $i(fccr, qr, k - qar) 
165 N ^ iuj~ E qT - E k ^ q — 

qr H H 

h 2 c \ 2 2 y-^^l(ka,qT,k + qaf) 
~~L6SN ^ iuj + E qT + E k+qW r 



, (B.3) 



, (B.4) 



E k a = E ka + Re S 1/S (fc, E ka + tO) . 



(B.l) 



where a = — a denotes a sign change such that ar 
—err, and the functions $1 and $2 are defined as 



§\{kxCri,k 2 (T 2 ,k 3 <J 3 ) = <Ti7i(ui ffl + aiVl<r t )(a 3 U2<r 2 V3o- 3 + a- 2 u 3a3 v 2 a 2 ) 

+ 0"272 (w 2c r 2 +C r 2«2 C r 2 ) ("10-1^3,73 + OzaiV 3 a 3 V lai ) 

+ 0-373 (u 3<73 + <T 3 y ao . 3 )(u 2o . 2 u lo . 1 + a 1 a 2 VxaiV 2 „ 2 ) , 

§ 2 (k 1 a 1 ,k 2 a 2 ,k 3 (T 3 ) = (J^xiu^, + a 1 v lai )(a 2 u 3a3 v 2rT2 + a 3 u 2a2 v 3a3 ) 

+ o- 2 7 2 (w2 C r 2 + fy 2 v 2<J2 ){a 1 u 3a3 v lcJl + a 3 u lai v 3a . 3 ) 

+ O-373 (U 3a . 3 + (7 3 V 3<J3 ) (aiU 2a2 V lcri + (T 2 U lai V 2 a 2 ) ■ 



(B.5) 



(B, 



The frequency independent 1/5-contributions to the self-energy are 
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-,1/S 
^3er 



-,1/S 

M<7 



(fe) 



(fe) 



25' 

25' 
ft, 



(u feo . + V ka )[— KiCt + K 2 n 2 — n 3 m 2 + CT7fc(— K 3 a + K 4 n 2 /2 — Kim 2 
cTUko-Vko-[K2m 2 - K 3 n 2 + 2a"/ k (n 2 a - Kin 2 + K i m 2 /2)} , 



c ™2 



m (« 2 - «1 + K3)[(Wfea + V ka )(l - <7J k ) - ZlkUko-Vko 



(B.7) 
(B.8) 



r 



la 



0.98 



0.96 



0.94 




FIG. 4: Evolution of the spin-wave velocity of the gap- 
less magnon as a function of the external magnetic held for 
S — 1/2. The full line shows the spin-wave velocity obtained 
numerically from Eq. ()B.1|) normalized by the zero-field value 
Co ~ L16co for 5=1/2 (see Ref. Q)). The dashed line shows 
the prediction of Eq. (|5,16ap . Good agreement is obtained in 
the limit of vanishing fields which confirms that the leading 
field dependence is described by Eq. (|5.16a[) . 



While the self-energy (|B.2|) can be easily evaluated nu- 
merically, it is not very accessible for analytical treat- 
ments and the leading small field behavior of the spin- 
wave dispersion is not easily extracted from it. The 
equivalent expression Eq. (|4.7[) in the Hermitian field 
parametrization is more amenable to an analytical inves- 
tigation of the long wavelength physics. To calculate the 
self-energy given in Eq. (|B.2[) we performed a two dimen- 
sional integration and used an analytical continuation to 
real frequencies. Performing a numerical derivative with 
respect to the momentum k at the point in the Bril- 
louin zone where the dispersion is gapless finally yields 
the spin- wave velocity. In Fig. 0] we compare the numer- 
ically obtained spin- wave velocity of the gapless mode at 
small fields with the prediction of Eq. (|5.16a|) . At very 
small fields, the numerical solution indeed confirms the 
behavior given in Eq. (|5.16a|) . For slightly larger fields, 
corrections beyond the linear dependence are also visible. 
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